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Abstract 

We have made the first direct measurement of the parity- violating coupling of the Z° 
boson to the strange quark, A s , using ~550,000 e + e _ — ► Z°^hadrons events produced 
with a polarized electron beam and recorded by the SLD experiment. Z° — > ss events 
were tagged by the absence of B or D hadrons and the presence in each hemisphere of a 
high-momentum K ± or K®. From the polar angle distributions of the strangeness-signed 
thrust axis, we obtained A s = 0.895±0.066(s£ai.)±0.062(s?/st). The analyzing power and 
uu + dd background were constrained using the data. 
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The extent of parity violation in the electroweak coupling of the Z° boson to an 
elementary fermion / can be specified by the parameter Af = 2v faf/(Vf+aV), where Vf 
(a/) is the vector (axial- vector) Zff coupling. In the Standard Model (SM), universal 
couplings are expected for the leptons (A e = = A T ), the down-type quarks (Ad = 
A s = Ab) and the up-type quarks (A u = A c = A t ). Precise measurements of the Af 
allow stringent tests of the SM and sensitivity through radiative corrections to e.g.: the 
top quark and Higgs boson masses (A e ^ T ); new physics that affects primarily the right- 
handed couplings (A djS ^ b ); and new physics that couples more strongly to heavier quarks 
(deviations from universality). 

All except A t can be measured in e + e~ annihilations at the Z° resonance via forward- 
backward production asymmetries in Of, the polar angle of the outgoing / with respect 
to the incoming e~ beam. At the SLC, the e~ beam has longitudinal polarization P e , 
the e + beam is unpolarized, and the Born-level differential cross-section for the process 
e+e-^Z°^ff is: 

(l-A e P e )(l+x 2 )+2A f (A e -P e )x, (1) 

where the last term is antisymmetric in x = cosOf. Using both left- (P e < 0) and right- 
polarized (P e > 0) beams of magnitude \P e \, one can measure both the initial- (A e ) and 
final-state (Af) couplings directly fl], |j; f° r Pe = one can measure only their product, 
or A f FB = 3A e A f /4:. 

The most precisely measured coupling is A e , with a relative error of 1.3% 0, U, an d 
lepton universality is verified at the 8% level || . In the quark sector, several measurements 
of A), and A c that use properties of the leading B and D hadrons can be combined to 
yield precisions of 2.0% and 4.4%, respectively ||. However there are few measurements 
of A u , A d or A s H, |5[] because the leading particles in u, d and s jets are more difficult 
to identify experimentally; they have relatively low energy, are not unique to events of a 
particular flavor, and nonleading particles of the same species are produced in hadronic 
jets of all flavors. Furthermore, these aspects of jet fragmentation are not well measured, 
and previous indirect measurements either relied on imprecise constraints from their data 
(OPAL: A£ B = 0.040±0.073; Aj? B = 0.068±0.037 J§) or are model-dependent (DELPHI: 
AJ. B = 0.101±0.012 §). 

In this Letter we present the first direct measurement of A s . We used high-momentum 
K and K® to tag Z ^ ss events, and the charge to separate s jets from s jets. The 
heavy flavor (cc + bb) event background was suppressed by identifying B and D decay 
vertices. The uu + dd background was suppressed and the s-s separation enhanced by 
requiring an s/s-tag in each event hemisphere, reducing any model dependence. The 
remaining uu+dd background and the s-s separation were constrained using related ob- 
servables in the data. 

We used the sample of approximately 550,000 hadronic Z° decays recorded by the 
SLD § experiment at the SLAC Linear Collider, with (\P e \) = 0.735±0.005 §, from 1993- 
1998. Charged tracks were measured in the Central Drift Chamber (CDC) and the 
original (upgraded) Vertex Detector (VXD) [|[ in 26.5% (73.5%) of the data; the resolution 
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on the impact parameter d in the plane perpendicular to the beam direction, including 
the uncertainty on the interaction point, was =ll©70/(josin 3/2 fl) fim (8©29/(p sin ' U 
/^m), where p is the track momentum in GeV/c and 9 its polar angle with respect to the 
beamline. Tracks were identified as n^, or p/p in the Cherenkov Ring Imaging 
Detector (CRID) [0], which allowed the identification with high efficiency and purity of 
vr ± with 0.3 < p < 35 GeV/c, with 0.75 < p < 6 GeV/c or 9 < p < 35 GeV/c, and 
p/p with 0.75 <p< 6 GeV/c or 10<p<46 GeV/c JTof. The event thrust axis was 
calculated using energy clusters measured in the Liquid Argon Calorimeter . 

After selecting hadronic Z° decays |R|, we removed cc and 66 events by requiring 



no more than one well-measured |13j track with d/ad > 2.5 in the event. The efficiency 
for selecting light-flavor events with | cos 9 thrust I <0.71 and the VXD, CDC and CRID 
operational was estimated to be over 95%; the selected sample comprised 205,708 events, 
with an estimated contribution of 14.2% (3.4%) from cc (66) events. Such performance 
parameters were estimated from a detailed Monte Carlo (MC) simulation [|T^, |TJ] of the 
SLD based on the JETSET 7.4 JL5| event generator, tuned to reproduce many measured 



properties of hadronic Z° decays, including the momentum-dependent production of K ± , 
K°, K* and mesons. 

Each selected event was divided into two hemispheres by the plane perpendicular to 
the thrust axis, and in each hemisphere we searched for high-momentum strange particles 
K , and A°/A°. Candidate K ± tracks were required to have p>9 GeV/c, d<l mm, 
to extrapolate through an active region of the CRID gaseous radiator system, and to have 
a log-likelihood |L0| for the K ± hypothesis Ck that exceeded both C w and C p by at least 
3 units. For p>9 GeV/c, the estimated K ± selection efficiency (purity) was 48% (91.5%). 

Candidate K® ^7r + 7r~ and A°/A°^p7r~/p7r + decays were reconstructed as described 
in [ID, |T6j from tracks not identified as K^. We required p>5 GeV/c and a reconstructed 
invariant mass m nn or m p7r within two standard deviations of the K® or A mass. If CRID 
information was available for the p/p track in a A /A candidate, we required C p > C n ; 
otherwise we required that the A°/A° not be a K® candidate and that the flight distance 
exceed 10 times its uncertainty. The estimated A°/A° reconstruction efficiency (purity) 
was 12% (90.7%). These A°/A° were removed from the K® sample, for an estimated K® 
efficiency (purity) of 24% (90.7%). 

We considered only the selected strange particle with the highest momentum in each 
hemisphere (5.5% of those tagged contained more than one), and tagged the event as ss 
if one hemisphere contained a and the other contained either an oppositely charged 
K ± or a K^. The A°/A° tags provided a useful veto in multiply tagged hemispheres and 
important checks of the simulation; however their inclusion did not improve the total error 
on A s . The thrust axis, signed so as to point into the hemisphere containing (opposite) 
the K~ (K + ), was used as an estimate of the initial s-quark direction. Table [I] shows 
the number of events tagged in each mode, along with the predictions of the simulation, 
which are consistent. Also shown are the simulated ss event purities and analyzing powers 
a s = (N r — N w ) I (N r + N w ) , where N r (N w ) is the number of events in which the signed 
thrust axis pointed into the true s (s) hemisphere. 



3 



Mode 


# Events 
in Data 


MC 
Prediction 


ss 
Purity 


Analyzing 
Power 


K + K- 


1290 


1312 


0.73 


0.95 


K±K° S 


1580 


1617 


0.60 


0.70 


K+A°, K-A u 


219 


213 


0.66 


0.89 


A°A° 


17 


14 


0.57 


0.70 


A K° S ,A K° S 


193 


194 


0.50 


0.32 



Table 1: Summary of the selected event sample for the two tagging modes and the three 
cross-check modes. 



Figure [I] shows the distributions of the measured s-quark polar angle 8 S for the K + K~ 
and K ± Kg modes. In each case, production asymmetries of opposite sign and different 
magnitude for left- and right-polarized e~ beams are visible. The content of the largest 
| cos 9 S | bins is reduced by the detector acceptance. The estimated backgrounds (discussed 
below) are indicated: those from cc+bb events exhibit asymmetries of the same sign and 
similar magnitude to those of the signal, so the measured A s is largely insensitive to them; 
those from uu + dd events exhibit asymmetries of opposite sign, and A s is more sensitive 
to the associated uncertainties. 

A simultaneous maximum likelihood fit to these four distributions was performed using 
the function: 

L= if E N q {(l-A e P e )(l+x 2 k )+2(A e -P e )(l+5)a q A q x k }. (2) 

k=l q=udscb 

Here, the number of tagged qq events N q = N events R q e q , R q = T(Z°— >gg)/r(Z°^hadrons), 
e q is the tagging efficiency, a q is the analyzing power for tagging the q direction, and the 
correction for hard gluon radiation 5 = —0.013 was derived [T7| as in [TJJ. The values of 



the e q and a q depend on the tagging mode. World average values || of A e , A c , A b , R c 
and Rb were used, along with SM values of A u , Ad, R u , Rd and R s . Simulated values of 
e c , €b, a c and a& were used, as they depend primarily on measured quantities with well 
defined uncertainties. 

For the light flavors, the relevant parameter values were derived where possible from 
the data. The number of events N u +N d +N s was determined by subtracting the simulated 
N c and A^ from the total observed. The values of a s and the ratio (N u + No) /N s were 
constrained (see below) using the data; since the simulation was consistent with the data, 
the simulated values of a s were used and the simulated e u , and e s were scaled by a 
common factor to give the measured N u +N d +N s . The average a ud = (N u a u +N d ad) / (N u + 
Nd) can also be constrained from the data; however our constraint is less precise than 
the range — a s < a u , a d < 0, obtained by noting that a u (d) jet can produce a leading 
K + (K*° — > K + 7i~), giving a u (a,d)<0, but with an associated K~ or K°, and a K~ can 
be selected with reduced probability, giving |o u |(|ad|)< \a s \. We scaled the simulated 
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Figure 1: Measured s-quark polar angle distributions (dots) for selected events in the a,b) 
K + K~ and c,d) K ± K® modes, produced with a,c) left- and b,d) right-polarized electron 
beams. The histograms represent the result of a simultaneous fit to the four distributions, 
and the upper (lower) hatched areas indicate the estimated uu+dd (cc+bb) backgrounds. 



a u and by a common factor such that a u d = — a s /2 for each mode. The fit yielded 
t4 s = 0.895±0.066 (stat.). Histograms corresponding to this value are shown in fig. [I] and 
are consistent with the data; the binned \ 2 is 42 for 48 bins. 

We considered several sources of systematic uncertainty, summarized in table |[ The 
values of R c , R b , A c and A b were varied by the uncertainties on their world averages ||. 
A large number of quantities in the simulation of heavy flavor events and detector per- 
formance were varied as in |14] with negligible effect on the measured value of A s . The 
yield and analyzing power of true K ± from D (B) decays have been derived from SLD 
data in the context of a measurement ]18[ of A c (A b ), and our simulation reproduces them 
within the measurement errors. We applied corresponding relative variations of ±9% to 
e c , ±3.3% to e b , ±5% (15%) to a c and ±3.6% (4.4%) to a b for the K + K~ (K ± K®) mode. 
The sum in quadrature of the uncertainties due to heavy flavor background is listed in 
Table |2|; the largest contribution is from Sa c . 

The key to this measurement is the understanding of the light-flavor parameters, for 
which there are few experimental constraints, and these gave rise to the dominant sys- 
tematic uncertainties (Tj|. In order to minimize model dependence, we used our data 
to constrain the largest uncertainties in these parameters within the context of our sim- 
ulation, which reproduces existing measurements of relevant quantities such as leading 
particle production and strange-antistrange correlations [fTE] . 

To constrain the analyzing power a s , we note that to mistag an s jet as an s jet 
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we must either identify a true K + or misidentify a f + or p as a K + . A true high- 
momentum K + in an s jet must be produced in association with an antistrange particle, 
yielding a jet with three high-momentum particles of nonzero strangeness. In our data we 
found 61 hemispheres containing three selected and/or K®] the MC prediction of 67 is 
consistent. We quantified this as a constraint on a s by subtracting the simulated (c+c+6+6) 
contribution of 9.3, scaling by the simulated ratio (s+s) / (u+u+d+d+s+s) = 0.74, and 
comparing with the MC prediction for (s + s). Propagating the data and MC statistical 
errors yielded an 18.5% relative constraint on the wrong-sign fraction, w s = (1 — a s )/2, 
in s/s hemispheres. This constraint is not entirely model- independent but any further 
uncertainties are small compared with 18.5%. Assuming equal production of charged and 
neutral kaons, this procedure delivers a calibration of w s for both tagging modes, which we 
varied simultaneously by ±18.5% relative. To account for misidentified particles we varied 
the production of >9 GeV/c p and 7r + in s jets by ±100%, and varied the misidentification 
probability by its measured relative uncertainty of ±25% [TO]. The sum in quadrature of 
these three effects is shown in table |] and is dominated by the 3-kaon calibration. 

The relative uu + dd background B u d = (N u + Nd)/N s was constrained in a similar 
manner, by exploiting the fact that an even number of particles with nonzero strangeness 
must be produced in a u or d jet. The three quantities, the number Aq = 1262 of hemi- 
spheres in the data containing an identified K + K~ pair, A 2 = 983 hemispheres containing 
a K^K® pair, and A3 = 503 events with an identified K of the same charge in both 
hemispheres, constrain B u d in complementary ways: Aq and A2 are primarily sensitive 
to KK production in u/d jets; (Aq — N 2 ) to (ft production in s jets; and A 3 to these and 
also the production and misidentification of tt^ and p/p. Furthermore, all are sensitive to 
deviations from the assumed values of R u , Rd and R s . The MC predictions of Aq = 1218, 
A 2 = 1002 and A 3 = 559 are consistent, and relative constraints on B U( i of 4.6%, 5.1% 
and 8.1%, respectively, were derived. Since A 3 constrains the sum of all contributions, we 
varied B u d by ±8.1%. 

These quantities are also sensitive to a u and ad, however our limited event sample 
did not allow us to obtain a useful constraint. We therefore took — a s < a u ,a d < as 
hard limits and scaled a u and simultaneously such that a ud = — a s /2±a s /\/l2- This 
yielded the dominant systematic error on A s and is a quantity that must be understood 
experimentally before a more precise measurement can be made. Since the product A q a q 
appears in Eqn. |2], this is equivalent to varying A u and Ad down to half of their SM 
values and up to well over unity; we considered no additional variation of A u or Ad- The 
uncertainties listed in table § were added in quadrature to yield a total relative systematic 
error of ±0.069. 

Several systematic checks were also performed. Ad hoc corrections ]TJ| to the simula- 
tion of the kaon momentum distributions and identification efficiencies, and the charged 
track reconstruction efficiency and impact parameter resolution were removed and the 
analysis repeated; changes in the measured value of A s were much smaller than the sys- 
tematic error. We fitted each tagging mode separately, including those involving A tags, 
with consistent results. We repeated the analysis using all K , and all A°/A°, hemi- 
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Source 


Systematic 
variation 


SA S /A S 


TTeavv flavor hark ground 


See text 


0.014 


Correction for gluon radiation 


-0.013±0.006 


0.006 




u. i oojzu.uuo 


u.uuu 


TVTO statistics 

1V1 V_ / U CI U XkJ U 1 V. k .J 




0.014 


a s K + K~ 

K±K° S 


0.949±0.012 
0.701±0.060 


0.032 


(N u + N d )/N s K + K- 
(incl. (R u + R d )/R s ) K^K* 


0.190±0.015 
0.316±0.026 


0.021 


Q"U) O'd 

A u , Ad 


-a s /2±a s / v / l2 


0.054 


Total 




0.069 



Table 2: Summary of the systematic uncertainties. 



spheres with no tag required in the opposite hemisphere; results were consistent. This 
K ± analysis is similar to that in ||; it has a relative statistical precision of 0.03, but of 
0.18 systematic. 

In conclusion, we have made the first direct measurement of the parity-violating cou- 
pling of the Z° boson to the strange quark, 

A s = 0.895±0.066(stat.) ± 0. 062 {syst.), 

using high-momentum identified and K® to tag Z°^ss decays and determine the s- 
quark direction. Our high K identification efficiency allowed the use of a relatively high- 
purity, double-tagged event sample, and the extraction from the data of constraints on the 
analyzing power of the method and the uu+dd background, using events with same-charge 
double tags and jets with two or three identified kaons. This result is consistent with the 
Standard Model expectation, A s = 0.935, and with less precise, previous measurements 
of A S FB ||, |5|. It is also consistent with a recent world average 6-quark asymmetry, 
v4 6 = 0.881±0.018 0, providing a 10% test of down-type quark universality. 
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